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Abstract 

We show that a single-mode squeeze operator S(z) being an unitary operator with a 
purely continuous spectrum gives rise to a family of discrete real generalized eigenvalues. 
These eigenvalues are closely related to the spectral properties of S(z) and the corre- 
sponding generalized eigenvectors may be interpreted as resonant states well known in 
the scattering theory. It turns out that these states entirely characterize the action of 
S(z). This result is then generalized to N-mode squeezing. 

1 Introduction 

Squeezed states play a prominent role in the modern quantum optics, see e.g. [1, 2, 3]. They 
are quantum states for which no classical analog exists. Recently, they have drown a great 
deal of interest in connection with quantum teleportation. Squeezed states of light were 
successfully teleported in the experiment reported by the group of Furusawa [4] . Mathemat- 
ical properties of these states were investigated in a series of papers [5, 6, 7, 8, 9]. In the 
present paper we analyze the spectral properties of single- mode S(z) and iV-mode S^(Z) 
squeeze operators. Clearly, squeeze operators are unitary and hence their spectra define a 
subset of complex numbers with modulus one. It is easy to show that the spectrum of the 
squeeze operator operator is purely continuous and cover the entire unit circle on the com- 
plex plane. However, it is not the whole story. Actually, it is easy to note that S(z) displays 
two families of discrete real eigenvalues. Clearly, these eigenvalues are not proper, that is, 
the corresponding eigenvectors do not belong to the corresponding Hilbert space of square 
integrable functions. One of these families was reported in a series of papers by Jannussis et 
al. [10, 11]. However, the interpretation of this result was not clear (it was criticized by Ma 
et al. [12] who stressed that proper eigenvalues of the squeeze operator do not exist). In the 
present paper we are going to clarify this problem. In particular we show that S(z) being 
an unitary operator with a purely continuous spectrum does indeed possesses discrete eigen- 
values which are closely related to its spectral properties. The corresponding eigenvectors 
may be interpreted as resonant states well known in the scattering theory [13, 14]. We show 
that restricting to a suitable class of states (e.g. coherent states do belong to this class) the 
action of S(z) may be entirely characterized in terms of these discrete eigenvalues and the 
corresponding eigenvectors. This observation is then generalized to two-mode squeezing and 
finally to iV-mode squeezing. 
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2 Single-mode squeezing 

A single-mode squeeze operator is defined by 

S(z) = exp Q[^o t2 -z*a 2 ^\ > (2- 1 ) 

where z is a complex number and a (a' ) is the photon annihilation (creation) operator which 
obey the standard commutation relation [a, a 1 "] = 1. Clearly, S(z) may be represented as 
S(z) = exp(iH(z)), with 

H{z) = 1 (zat 2 - z*a 2 ) . (2.2) 

Now, to investigate spectral properties of H(z) let us consider a unitarily equivalent operator 
R^((f)H(z)R((f), where R(ip) is a single-mode rotation [15, 16, 17] 

R(<p) =exp(i(^a^a) . (2.3) 

One has 

R\<p)H{z)R{v)=H{ze- 2i v) , (2.4) 

and hence, for tp = 9/2, where z = re 10 , it shows that H(z) is unitarily equivalent to H(r). 
Introducing two quadratures x and p via 

x + ip t x-ip 
a = ^ r , a^-^, (2.5) 

one finds the following formula for H{r): 

H(r) = ~ (xp + px) . (2.6) 

Spectral properties of H(r) were recently investigated in [18] (see also [19]) in connection with 
quantum dissipation. Note, that the classical Hamilton equations implied by the Hamiltonian 
H = —rxp: 

x = —rx , p = rp , (2-7) 

describe the damping of x and pumping of p. This is a classical picture of the squeezing 
process. It turns out that H(r) has purely continuous spectrum covering the whole real line. 
Hence, it is clear that it does not have any proper eigenvalue. Using standard Schrodinger 
representation for x and p = —id/dx, the corresponding eigen-problem H(r)ip = Eip may be 
rewritten as follows 

x^(x) = - (ij + ^ 1>(x) . (2.8) 

Note, that H(r) is parity invariant and hence each generalized eigenvalue E G R is doubly 
degenerated. Therefore, two independent solutions of (2.8) are given by 

^(x) = ^=x-W^ , (2.9) 
V 27rr 
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where x± are distributions denned as follows [20] (see also [21]): 

4-(f 111 . ^ ■■= \ ,i III . (-0, 



^ x < ' - " \ |x| A x < ' 

with A € C (basic properties of x± are collected in [18]). These generalized eigenvectors ip E 
are complete 

Ji^)^(x')dE = S(x -x 1 ) , (2.11) 

and (5-normalized 

J Vf(x)Vf 2 (x) dx = ^ - £ 2 ) . (2.12) 
Hence they give rise to the following spectral resolution of H(r): 

ff(r) = W El^Xi^ldE, (2.13) 

and the corresponding spectral resolutions of squeeze operator S(r) immediately follows 

S(r) = J2 [ e iE \^ E )(^ E \dE. (2.14) 

Now, let us observe that FH = —HF, where F denotes the Fourier transformation. Hence, 
if H(r)ip E = Ei[) E , then H(r)F[^ E ] = EF^~ E ]. Therefore, the family F[ipg E ] defines 
another system of complete and 5-normalized generalized eigenvectors of H(r). Note that 
the action of S(r) is defined by 

S(r)V>(x) = e" r /V(e" r x) , (2.15) 

and its Fourier transform 

F[S{r)4>](p) = e r/2 F[i;}(e r p) , (2.16) 

due to FS(r) = S(—r)F, that is, if the fluctuations of p are reduced then the fluctuation of 
x are amplified and vice versa. 

3 Discrete real eigenvalues of S(r) 

Surprisingly, apart from the continuous spectrum H(r) gives rise to the following families of 
complex discrete eigenvalues [18] 

H{r)ft = ±E n ft , (3.1) 

where 

E n = ir(n+^\ , (3.2) 
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and 

(_-\\n r n 

f-(x) = 1^ «W(x) , /+(*) = -= . (3.3) 
VI! V^! 

Interestingly they satisfy the following properties: 

/oo 
/n 0*0 fmi x ) dx = 5 nm , (3.4) 
-oo 

and 

oo 

Y.fn{x)f-{x') = 5{x-x') . (3.5) 

n=0 

It implies that 

S(r)ft = e ±iE - ft , (3.6) 
which shows that S(r) displays two families of purely real generalized eigenvalues 

r " 



st = ex P 



±r ( n + ^ 



(3.7) 



A family s+ was already derived by Jannussis et al., see e.g. formula (2.3) in [11], but they 
overlooked the second one s~. How to interpret these eigenvalues? It turns out that one 
recover E n and by studying a continuation of generalized eigenvectors ip± and F[\j}± E ] 
into the energy complex plane E G C. Both -0 s and F[ip~ E ] display singular behavior when E 
is complex: ip± has simple poles at E = —E n , whereas F[ip± E ] has simple poles at E = +E n , 
with E n defined in (3.2). Moreover, their residues correspond, up to numerical factors, to the 
eigenvectors f^: 

Res(Vf(x);-£„) ~ /" , (3.8) 

and 

Res(F[^ E (x)];+E n ) ~ /+ . (3.9) 

Such eigenvectors are well known in scattering theory as resonant states, see e.g. [13] and 
references therein. In the so called rigged Hilbert space to quantum mechanics these states 
are also called Gamov vectors [14]. To see the connection with the scattering theory let us 
observe that under the following canonical transformation: 

rQ-P rQ + P . 

X = ^r> p= ^r> (3 - 10) 



H(r) transforms into the unitarily equivalent operator [22] 

H(r) — H io = ^(P 2 -r 2 Q 2 ) , (3.11) 

which represents the Hamiltonian of the called an inverted or reversed oscillator (or equiva- 
lently a potential barrier '— r 2 Q 2 /2') and it was studied by several authors in various contexts 
[23, 24, 25, 26, 27, 28]. An inverted oscillator H m corresponds to the harmonic oscillator with a 
purely imaginary frequency u> = ±ir and hence the harmonic oscillator spectrum l u(n+ 1/2)' 
implies l ±ir(n + 1/2)' as generalized eigenvalues of H w . 
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4 A new representation of S(r) 

Interestingly, the action of S(r) may be entirely characterized in terms of and s^- Indeed, 
consider a space T> of smooth functions ip = ip(x) with compact supports, i.e. ip{x) = for 
\x\ > a for some positive a (depending upon chosen -ip), see e.g. [29]. Clearly, V defines a 
subspace of square integrable functions L 2 (R). Moreover, let Z = F[V], that is, ip € Z if 
tp = F[<fi] for some 4> £ V. It turns out [29] that V and Z are isomorphic and V n Z = 0. 
Now, any function 4> from Z may be expanded into Taylor series and hence 

oo i(n) (n\ 00 
n=0 ' n=0 

On the other hand, for any <j) G T>, its Fourier transform F[<f>] 6 Z, and 

0(x) = 1 ( e ik *F[<t>]{k)dk = -L fe«^fffS F& 

oo oo 

= ^F[/ n + ]W(/;W]) = ^/ n -W(/;w. (4.2) 

n=0 n=0 

Hence, we have two decompositions of the identity operator 

oo 

1 = E l/n)</n I °n Z, (4.3) 
n=0 

and 

oo 

l = El/rT)</^l °n P. (4.4) 

n=0 

It implies the following representations of the squeeze operator S(r): 

oo 

S(r) = J2 s n\fn)(fn\ on Z, (4.5) 

n=0 

and 

oo 

S(r)=X>+|/-></+| 011 P - ( 46 ) 

n=0 

It should be stressed that the above formulae for S(r) are not spectral decompositions and 
they valid only on Z and V, respectively (its spectral decomposition is given in formula 
(2.14)). Note, that S(r) maps Z into V and using (4.5) one has 

00 

S\r) = Y J S-\f-){ft\ = S{-r) on V . (4.7) 

n=0 
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Conversely, S(r) represented by (4.6) maps T> into Z and S'(r) = S{—r) on Z. It shows 
that squeezing of x (p) corresponds to amplifying of p (x). Clearly, a general quantum state 
ijj G L 2 (R) belongs neither to V nor to Z. An example of quantum states belonging to Z 
is a family of Glauber coherent states \a). Consider e.g. a squeezed vacuum S(r)ipo, where 
i/jq(x) = ir~ 1 / 4: e~ x I 2 . One has 

and 

S(r)Mx) = —rnr Y r = ftJ?) = —rn^ Y ^— ftS^x) . (4.9) 

The similar formulae hold for S(r)\a). 



5 Two-mode squeezing 

Consider now a two-mode squeeze operator [15, 16, 17] 

S2(z) = exp ^za{(4 — z*aia,2*j , (5-1) 

where a\ and are creation and annihilation operators for two modes k = 1,2. Introducing 
2-dimensional vectors 

a T = ( ai ,a 2 ), (a t ) T = (4,4) , (5.2) 

one finds 

S2{z) = exp {z{a}) T <j\g) — z*a r aiaj , (5-3) 
where o\ stands for the corresponding Pauli matrix. Using well-known relation 

e i7T/ia2 a 1 e- i7T/4a2 =a 3 , (5.4) 

one obtains 

£ 2 ( z) e -W4<x 2 = exp ^(at) T <73 a t - z*a T a 3 a) 

= S (1) (z)S (2) (-z) , (5.5) 

where S^* 1 denotes a single-mode (fc) squeeze operator. Now, since (z) is unitarily equiv- 
alent to S( k \r) a two-mode squeeze operator S^z) is unitarily equivalent to S^^r)^ 2 ^— r). 
Hence, the spectral properties of S2 {z) easily follows. In particular defining the space V 
of smooth functions tp = 4 } { x i, x 2) with compact supports and Z = F[V] one obtains the 
following representations: 



00 



S (i) {r)S (2) { _ r)= £ s - m |/+ m )(/-J on Z, (5.6) 



nm=0 
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and 



S (i) {r)S (2) { _ r)= J- s+ m \f- m )(f+J on V, (5.7) 

nm=0 

where 

f± m (x 1 ,x 2 ) = fn (zi)/™(x 2 ) , (5.8) 



and 

stm = e ±rin ~ m) • (5.9) 

Jannussis et al. [11] claimed that the eigenvalues of S 2 (z) are given by e 2 ( m - n ) ) se e e.g. 
formula (5.7) in [11]. Their result has the similar form as s~ m but of course it is incorrect. 
Note that eigenvalues of [11] do not depend upon the squeezing parameter z as was already 
observed in [12]. 

6 N-mode squeezing 

Following [30] one defines an TV-mode squeeze operator 

S N (Z) = exp Q (g?) T Za ] - Z f a) (6.1) 

where Z is an TV x TV symmetric (complex) matrix and 

a T = (a±,a 2 , ■ ■ -,a N ) . 
Defining an TV-mode rotation operator 

R N {§) = exp (i(a f ) T 8a) , (6.2) 

with being an TV x TV hermitian matrix, one shows [17, 30] 

4,(8) S N (Z) R N @) = S N (e~ S Z e'^ T ) . (6.3) 

Now, by a suitable choice of $> one obtains 

e~ S Z = Z D , (6.4) 

where Zd is a diagonal matrix, i.e. (^d)a;/ = Zk$ki- Hence, an TV-mode squeeze operator 
Sn(Z) is unitarily equivalent to 

4(8) S N (Z) R N ($) = SM( Z1 ) sW(z 2 ) ■ ■ ■ S (N \z N ) , (6.5) 

and therefore its properties are entirely governed by the properties of the single-mode squeeze 
operator S(z). In particular Sn(Z) gives rise to a discrete family of generalized eigenvalues 
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being combinations of generalized eigenvalues of S^ k \z k ). Defining the corresponding sub- 
spaces T> and Z in the Hilbert space L 2 (W N ) one easily finds 

oo 

RU$)S N (Z)R N ($)= S n 1 ...n N \ti 1 ...n N )(fn 1 ...n N \ «n Z, (6.6) 

ni,...,njv=0 

and 

oo 

RU*)Sn(Z)R N {*)= S n 1 ...n N \fn 1 ...n N )(ft..n N \ «n V, (6.7) 

m,...,n N =o 

where 

fn 1 ...n N ( x U---' X N)= fn^Xl) ■ ■ ■ fn N ( x N) , (6.8) 

and 

»ni...njv = exp j± 

with r k = \z k \. 
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